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ABSTRACT 


The most significant limitation of the imaging SAR in rain measurements is the ground 
return coupled to the rain cell. Here we report a study of the possibility of using the X-SAR and 
the C-band channel of SIR-C for rain measurement. Earlier signal-to-clutter calculations rule out 
the use of X-SAR at steeper off-vertical angles of incidence ( i.e., 20 < 0 < 50). Only rain rates 
> 30 mm/hr at angles of incidence > 60° showed good signal-to-clutter ratio (SCR). 

This study involved calculations at vertical incidence. There is adequate signal-to-noise 
ratio (SNR) at vertical incidence, but the presence of high-range side-lobe levels leads to small 
SCR for measurement over oceans at both X and C bands. For larger rain thickness (> 2 km), 
the SCR gets better and smaller rain rates (> lOmm/hr) can be measured. However, rain 
measurements over forests seem to be feasible at nadir even for smaller rain thickness (< 2 km). 
We conclude that X band may be usable over the forest at vertical incidence to measure rain 
rates > 5 mm/hr even for shallow rain thickness and over ocean for large rain thickness. 


1.0 INTRODUCTION 

The detection and measurement of precipitation from space platforms is important in 
various meteorological applications. A precipitation detection and measurement experiment was 
accepted for the SIR-C/X-SAR mission. The proposed experiment will use the SIR-C /X-SAR 
sensors, to measure rain quantitatively for moderate to intense rain rates ( > 10 mm/hr). X-SAR 
achieves adequate SNR at all angles of incidence for all but the lightest rains. The C-band SIR- 
C gets reasonable SNR at all angles for moderate to heavy rains. However, the SCR is so poor 
at C band that it cannot be used over any reasonable ground surface. Earlier it was concluded 
that it would be possible to measure rain rates in excess of 30 mm/hr at angles of incidence 
greater than 60°. For this study we use the backscatter method for rain-rate measurement. The 
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backscatter algorithm is the simplest of algorithms; however its use without corrections for 
attenuation leads to considerable error in the rain rate estimation (Goldhirsh, 1988). 

2.0 SYNTHETIC-APERTURE RADAR RESOLUTION IN RAIN 

When observing a target, the maximum possible aperture that can be synthesized during 
time T a for which the target is illuminated is 

L s = UT a (1) 

where U is the platform velocity. Atlas and Moore (1987), showed that the expected Doppler 
spread from rain permits only a relatively short integration time. A Synthetic-aperture radar 
(SAR) has better resolution than a real-aperture radar only when 


where L n is the length of the along-track real aperture, U is the platform velocity and o vs is the 
standard deviation of the Doppler velocity spectrum. The critical value of o vs is 3.75 m/s for X 
band and 6.81 m/s for C band. The SAR can use the full synthetic-aperture only when the 
targets in the pulse volume are correlated (i.e., their phase positions remain fixed). This means 
that the effective beamwidth, with rain as the target, increases with the increase in o vs , the spread 
of the target Doppler spectrum. 

2.1 GEOMETRY 

The echo power for the SAR is identical to that for the conventional real aperture radar, 
using the same antenna. Hence, will use the same geometry as that of real-aperture for power 
calculations. Fig. 1 shows the geometry for the vertical incidence. The rain volume may be 
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Fig.l The geometry with the shuttle looking at nadir. 
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divided into three regions. 

Region 1. where there is partial beam filling due to the pulse intercepting the cloud top. 
Region 2. where the cells are assumed to be completely rain filled. 

Region 3. where the rain cells intercept the surface. 

Here we consider calculations only for those cells that are completely rain filled (ideal 
rain cells). The volume of these rain cells is 

v n = -f P* Pv r R R n 2 (3) 

where R n is the range to the cell in question, (3 h and p v are the horizontal and the vertical beam 
widths, and r R is the vertical resolution. The factor of tc/4 accounts for the elliptical shape of the 
beam. The ground area coupled to each cell except the first can be calculated, assuming the area 
is approximately rectangular (see Fig. 2), as 

A n = 2 P /3 R a ( x{n) - x(n-i) ) (4) 

and the ground area coupled to the first cell is given by 

= 2 P/i R n x(l) ( 5 ) 

where 

x(n) = R n sin (0 n ) (6) 

0„ = cos- 1 (A) (7) 

K n 

and 

R n = h s + n r R (8) 

We define h s as the shuttle altitude, R n as the range to the particular cell in question and n as the 
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cell number ranging from 1 to m. To calculate m, the total number of cells coupled to the 
ground, we use the relation 


R s - h s 
m = — 

1 R 

where R s is the maximum slant range to the ground (see Fig. 3) 


R 


S 


COS ( P v / 2 ) 


(9) 


( 10 ) 


3.0 MODELS USED 

For vertical incidence, the ocean return at X band is assumed to have a scattering 
coefficient o° of 12 dB (Schroeder et al., 1985) and forests to have a a° of -6 dB ( Ulaby et al., 
1986) . Further, due to the large magnitude of the returns, backscatter modification due to wind 
speeds and impinging rain drops is ignored. 

Radar meteorologists have derived an empirical power law, the Z-R relationship, from the 
available data and Mie scattering (Battan,1973) . That is 

z = a R*' I 11 ’ 

where a and fr, have been found to depend on the dropsize distribution, rain type, rain rate, 
frequency, and the region. Although it is clear that the Z-R relation is useful, the values of the 
coefficient must be selected judiciously. The coefficients differ markedly at small and large 
values of rain rate R,, hence a single Z-R relationship cannot be used. For modest and small 
rain rates, it is resonable to select a around 300 and b { around 1.5 (Atlas and Moore 1987). 
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F»g.3 The geometry to calculate the number of cells coupled with the surface echo. 
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4.0 SIGNAL-TO-NOISE RATIO 


From the basic radar equation 


P, G 2 X 2 o v 

p — C 1 

r (4 tc ) 3 i? 4 


( 12 ) 


where the variables have their usual meaning, we obtain an equation for SNR. For the vertical 
incidence angles, using (3) for volume of the ideal rain cells, (12) may be written as 


where 


SNR n = 


K T B F 


1024 n 2 R‘ 


(13) 


n 


tc 5 300 J? 1 - 5 ! k\ 1Q~ 10 
X 4 


(14) 


Here G is the gain, B is the bandwidth, L s is the system loss, 1^, is the signal loss due to 
absorption by the melting ice layer, T| is the rain reflectivity, F is the receiver noise figure, and 
the other parameters have their usual meaning. The extinction coefficient [see Appendix I] is 


where 


k er [ Rr } 



a[Rr } 


2 


(h d i 


cos (6) -1 k 
cos ( 8 ) 



a [J?r] 


b Rr c _ 

8.66 


(15) 


(16) 


The value of b and c depend on the frequency used (see Table 1), h d is the shuttle’s height above 
the top of the rain cell and h k is the distance of the rain cell from the rain top. Using (13) we 
plot the SNR in dB vs rain rate in Fig. 4. 

Only calculations for the cells that are entirely rain filled are presented here. For partially 
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Fig.4 Plot of the SNR in dB vs Rain rate in mm/hr for vertical incidence. 
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rain-filled cells, the SNR decreases in proportion to the rain-filled volume. 


5.0 SIGNAL-TO-CLUTTER RATIO. 

The backscatter method uses the returned echo from the rain to compute the rain rate. 
For vertical incidence the surface returns, coupled by the sidelobes, arrive simultaneously with 
the desired echo and increase in magnitude as we approach cells near the surface. For the 
backscatter method to work, the echo from the rain should be substantially above the surface 
echo to allow separating them. Hence, we compute the signal-to-clutter ratio. This is defined 
as the ratio of the desired signal (rain return) to the surface return. Since the surface and the rain 
are excited by the same transmitter, one need not use the radar equation in its entirety. 

The scattering coefficient of the rain, when X is in centimeters, is 


a= — Z \k\ 2 1CT 10 
r A 4 


(17) 


Using the appropriate volume from (3) and the extinction coefficient from (15) 


o.. - a. vr 


V I 


Pv 


(18) 


and the appropriate area coupled to the rain cell from (4) in 


o = o 0 (A n ) , 


(19) 


where the magnitude of ct° depends on the type of surface. We compute the SCR 


Ao rf „ = a v - 

US V dB ffdB 


( 20 ) 


where 


°v dB = 10 lo 9(° v ) 


( 21 ) 


and 
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a ffdB = 10 1 og(a ffT ) 


( 22 ) 


For the ideal case, the clutter coupled to the rain returns at vertical incidence is negligible. 
However in practice the surface return from range sidelobes arrives at the receiver at the same 
instant as the desired rain echo. The sidelobe levels with a raised cos 2 window are low, with the 
first side lobe about 21 dB down. Due to the very large magnitude of the surface clutter and due 
to the fact that the side lobes fall too slowly (see Fig. 5), almost all rain cells have significant 
clutter coupled to them. From Fig. 5 it is seen that the magnitude of the sidelobes that are 
coupled to the surface echo depends on the cell for which the calculations are made. Further, 
even the sidelobes that are coupled to the ground return change in magnitude with the cell under 
consideration. To account for this, using (9) and the cell number under consideration, it is 
possible to obtain the exact number of side lobes coupled to the ground. Their magnitude Sl dB (n), 
can be obtained from the figure. The weighted clutter coupled to the cell under consideration 
is obtained as 

Q gr U2 r ) = s a g ( R r ,n ) 10 01 SJ -® (n) (23) 

n 

where 


o g (R r ,n) 


A n K gi [Rr] 10 

B„ 


0.1 o„ 


(24) 


Using (20) plots of the SCR for various rain thicknesses over ocean are compared with 
SCR over forest for rain thickness of 1 km at X band (see Fig. 6(a)). Fig. 6(b) shows the 
comparison between the SCRs obtained using a raised cos 2 window and an ideal cos 2 window 
over the ocean at X band. Fig, 6(c) and 6(d) show corresponding comparisons at C band. From 
the figures it is seen that the X-band and C-band SCR with a raised cos 2 window is poor 
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Tablc.l The attenuation coefficients at various frequencies may be picked 


off from here to be used in the relation discussed. 


(Lin, 1975; Morrison and Cross, 1974) 
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over ocean over forest 


Fig.6a. Plot of X -band SCR in dB vs Rain rate in mm/hr over ocean and 
over forest for the given rain thickness for raised cos 2 weighting. 
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Rain Rate in mm 



Fig.6b Comparison of the X-band SCR over ocean using raised cos 2 and 
ideal cos 2 weighting for the given rain thickness. 
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SCR in dB 



raised cos~2 ideal cos A 2 


Fig.6d Comparison of the C-band SCR over the ocean for a raised cos 2 
and a ideal cos 2 weighting for the given rain thickness. 
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over the ocean for shallow rain thickness. However over the forest the SCR appears to be 
reasonable even for shallow rain thickness with a raised cos 2 window. Further for both X and 
C bands the SCR is expected to be good even for shallow rain thickness over ocean if an ideal 
cos 2 weighting could be achieved. 

6.0 CONCLUSIONS 

Though the SAR has a good SNR at all angles of incidence, it has inadequate SCR under 
most conditions. For the backscatter algorithm to work, the return from rain must have a 
magnitude greater than the surface return. From the plots it is seen that the SCR (which gives 
a measure by how much the rain return lies above the ground return) at nadir is very low for 
small rain rates over oceans (see Fig. 6(a)). This implies that the backscatter method cannot be 
used over ocean for rain rates < 30 mm for shallow rain thickness using a raised cos 2 weighting. 
However, the SCR for forest at vertical incidence is adequate for rain rates > 5 mm/hr even for 
rain thickness of 1 km at X band and > 10 mm/hr at C band. For low rain thickness, 
measurements over the ocean are not feasible for both vertical incidence and off-vertical 
incidence except at very high rain rates. However, for rain thickness > 2 km, measurement over 
ocean at X band for rain rates > 10 mm is possible at nadir. For C band the SCR is very poor 
even at high rain rates as seen in Fig. 6(b). 

Since far-out sidelobes are the problem at vertical incidence, it is suggested that the 
sidelobe level be reduced by an additional 10 dB during processing so that the experiment may 
be conducted over the ocean at nadir for shallow rain thickness. Here the absolute level of the 
sidelobes is not the issue and any weighting function that gives sidelobes which rapidly decrease 
in magnitude is recommended. Clearly with ideal cos 2 weighing the SCR is expected to be very 
good even at shallow rain rates (Fig. 6(b), Fig. 6(d)). However, over forest X-band and C-band 
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vertical incidence results are expected to be good even with the present weighting. This does not 
mean that experiment is feasible over forest at all angles because the clutter return from forest 
is higher than that of the ocean at incidence angles far from the vertical. 

It is believed that the rain retrieval algorithm could work for low rain rates at shallow rain 
thickness if it is possible to have an independent estimate of the signal and clutter and the clutter 
alone. By subtracting the average clutter estimate from the average rain signal-plus-clutter 
estimate it is possible to implement the rain retrieval algorithm. This will be discussed in a 


separate report. 
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APPENDIX - I 



Fig.a The geometry of the pulse volume showing that the signal from different parts of the pulse 
volume passes through different rain thickness 
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From Fig. a ii is seen that the volume occupied by the rain above each rain cell changes from 
one end of the cell to the other. This is due to the fact that pulse volume is curved where as the 
rain top is not. Hence the signal from the cell edge suffers less attenuation than the center. To 
account for this it is necessary to integrate the attenuation over the entire pulse volume. Here 
it is necessary to make the assumption that the cell has constant attenuation in the azimuth 

direction. The rain returns involve 


Jv 

o i Pa r, I X* e '’ R ‘ dR (a) 

h 

2 

we may write R = (h + n r R ) where n is the cell number under consideration, r R is the range 
resolution and h is the height to the shuttle from the rain top. In order to account for the 
variation of the distance travelled by the signal in the rain volume it is assumed that the 
beamwidth can be divided into number of small angles. From Fig. a 


R 


Z 


h + n r R 


h 

COS0 


(b) 


where 0 varies from - pv/2 to pv/2. Substituting (b) in (a) for the rain returns we have 


o v -jP* r R ( h + n z r )2 f e 
4 K 


-a ( 


h (cos 8 - l) 
COS 0 ~~ ' 


n r R ) 


de 


(c) 


Here the attenuation coefficient 


k« = / 

P 


“ ' ~ <C cose ~ * n Zr) & 0 


(d) 
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